Development of microsatellite markers for sister species Linum suffruticosum and Linum tenuifolium in their overlapping ranges.

BACKGROUND: Microsatellite markers were developed for distylous Linum suffruticosum and tested in the monomorphic sister species Linum tenuifolium. These species are perennial herbs endemic to the western and northwestern Mediterranean, respectively, with a partially overlapping distribution area. METHODS AND RESULTS: We developed 12 microsatellite markers for L. suffruticosum using next generation sequencing, and assessed their polymorphism and genetic diversity in 152 individuals from seven natural populations. The markers displayed high polymorphism, with two to 16 alleles per locus and population, and average observed and expected heterozygosities of 0.833 and 0.692, respectively. All loci amplified successfully in the sister species L. tenuifolium, and 150 individuals from seven populations were also screened. The polymorphism exhibited was high, with two to ten alleles per locus and population, and average observed and expected heterozygosities of 0.77 and 0.62, respectively. CONCLUSIONS: The microsatellite markers identified in L. suffruticosum and tested in L. tenuifolium are a powerful tool to facilitate future investigations of the population genetics, mating patterns and hybridization between both Linum species in their contact zone.

Genomic analyses of the Linum distyly supergene reveal convergent evolution at the molecular level.

Supergenes govern multi-trait-balanced polymorphisms in a wide range of systems; however, our understanding of their origins and evolution remains incomplete. The reciprocal placement of stigmas and anthers in pin and thrum floral morphs of distylous species constitutes an iconic example of a balanced polymorphism governed by a supergene, the distyly S-locus. Recent studies have shown that the Primula and Turnera distyly supergenes are both hemizygous in thrums, but it remains unknown whether hemizygosity is pervasive among distyly S-loci. As hemizygosity has major consequences for supergene evolution and loss, clarifying whether this genetic architecture is shared among distylous species is critical. Here, we have characterized the genetic architecture and evolution of the distyly supergene in Linum by generating a chromosome-level genome assembly of Linum tenue, followed by the identification of the S-locus using population genomic data. We show that hemizygosity and thrum-specific expression of S-linked genes, including a pistil-expressed candidate gene for style length, are major features of the Linum S-locus. Structural variation is likely instrumental for recombination suppression, and although the non-recombining dominant haplotype has accumulated transposable elements, S-linked genes are not under relaxed purifying selection. Our findings reveal remarkable convergence in the genetic architecture and evolution of independently derived distyly supergenes, provide a counterexample to classic inversion-based supergenes, and shed new light on the origin and maintenance of an iconic floral polymorphism.

What can cold-induced transcriptomes of Arctic Brassicaceae tell us about the evolution of cold tolerance?

Little is known about the evolution of cold tolerance in polar plant species and how they differ from temperate relatives. To gain insight into their biology and the evolution of cold tolerance, we compared the molecular basis of cold response in three Arctic Brassicaceae species. We conducted a comparative time series experiment to examine transcriptional responses to low temperature. RNA was sampled at 22°C, and after 3, 6, and 24 at 2°C. We then identified sets of genes that were differentially expressed in response to cold and compared them between species, as well as to published data from the temperate Arabidopsis thaliana. Most differentially expressed genes were species-specific, but a significant portion of the cold response was also shared among species. Among thousands of differentially expressed genes, ~200 were shared among the three Arctic species and A. thaliana, while ~100 were exclusively shared among the three Arctic species. Our results show that cold response differs markedly between Arctic Brassicaceae species, but probably builds on a conserved basis found across the family. They also confirm that highly polygenic traits such as cold tolerance may show little repeatability in their patterns of adaptation.

Mutation accumulation opposes polymorphism: supergenes and the curious case of balanced lethals.

Supergenes offer spectacular examples of long-term balancing selection in nature, but their origin and maintenance remain a mystery. Reduced recombination between arrangements, a critical aspect of many supergenes, protects adaptive multi-trait phenotypes but can lead to mutation accumulation. Mutation accumulation can stabilize the system through the emergence of associative overdominance (AOD), destabilize the system, or lead to new evolutionary outcomes. One outcome is the formation of maladaptive balanced lethal systems, where only heterozygotes remain viable and reproduce. We investigated the conditions under which these different outcomes occur, assuming a scenario of introgression after divergence. We found that AOD aided the invasion of a new supergene arrangement and the establishment of a polymorphism. However, this polymorphism was easily destabilized by further mutation accumulation, which was often asymmetric, disrupting the quasi-equilibrium state. Mechanisms that accelerated degeneration tended to amplify asymmetric mutation accumulation between the supergene arrangements and vice-versa. As the evolution of balanced lethal systems requires symmetric degeneration of both arrangements, this leaves only restricted conditions for their evolution, namely small population sizes and low rates of gene conversion. The dichotomy between the persistence of polymorphism and degeneration of supergene arrangements likely underlies the rarity of balanced lethal systems in nature. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Genomic Signatures of Sexual Selection on Pollen-Expressed Genes in Arabis alpina.

Fertilization in angiosperms involves the germination of pollen on the stigma, followed by the extrusion of a pollen tube that elongates through the style and delivers two sperm cells to the embryo sac. Sexual selection could occur throughout this process when male gametophytes compete for fertilization. The strength of sexual selection during pollen competition should be affected by the number of genotypes deposited on the stigma. As increased self-fertilization reduces the number of mating partners, and the genetic diversity and heterozygosity of populations, it should thereby reduce the intensity of sexual selection during pollen competition. Despite the prevalence of mating system shifts, few studies have directly compared the molecular signatures of sexual selection during pollen competition in populations with different mating systems. Here we analyzed whole-genome sequences from natural populations of Arabis alpina, a species showing mating system variation across its distribution, to test whether shifts from cross- to self-fertilization result in molecular signatures consistent with sexual selection on genes involved in pollen competition. We found evidence for efficient purifying selection on genes expressed in vegetative pollen, and overall weaker selection on sperm-expressed genes. This pattern was robust when controlling for gene expression level and specificity. In agreement with the expectation that sexual selection intensifies under cross-fertilization, we found that the efficacy of purifying selection on male gametophyte-expressed genes was significantly stronger in genetically more diverse and outbred populations. Our results show that intra-sexual competition shapes the evolution of pollen-expressed genes, and that its strength fades with increasing self-fertilization rates.

Genome assemblies of three closely related leaf beetle species (Galerucella spp.).

Galerucella (Coleoptera: Chrysomelidae) is a leaf beetle genus that has been extensively used for ecological and evolutionary studies. It has also been used as biological control agent against invading purple loosestrife in North America, with large effects on biodiversity. Here, we report genome assembly and annotation of three closely related Galerucella species: G. calmariensis, G. pusilla, and G. tenella. The three assemblies have a genome size ranging from 460 to 588 Mbp, with N50 from 31,588 to 79,674 kbp, containing 29,202 to 40,929 scaffolds. Using an ab initio evidence-driven approach, 30,302 to 33,794 protein-coding genes were identified and functionally annotated. These draft genomes will contribute to the understanding of host-parasitoid interactions, evolutionary comparisons of leaf beetle species and future population genomics studies.

Selection on Accessible Chromatin Regions in Capsella grandiflora.

Accurate estimates of genome-wide rates and fitness effects of new mutations are essential for an improved understanding of molecular evolutionary processes. Although eukaryotic genomes generally contain a large noncoding fraction, functional noncoding regions and fitness effects of mutations in such regions are still incompletely characterized. A promising approach to characterize functional noncoding regions relies on identifying accessible chromatin regions (ACRs) tightly associated with regulatory DNA. Here, we applied this approach to identify and estimate selection on ACRs in Capsella grandiflora, a crucifer species ideal for population genomic quantification of selection due to its favorable population demography. We describe a population-wide ACR distribution based on ATAC-seq data for leaf samples of 16 individuals from a natural population. We use population genomic methods to estimate fitness effects and proportions of positively selected fixations (α) in ACRs and find that intergenic ACRs harbor a considerable fraction of weakly deleterious new mutations, as well as a significantly higher proportion of strongly deleterious mutations than comparable inaccessible intergenic regions. ACRs are enriched for expression quantitative trait loci (eQTL) and depleted of transposable element insertions, as expected if intergenic ACRs are under selection because they harbor regulatory regions. By integrating empirical identification of intergenic ACRs with analyses of eQTL and population genomic analyses of selection, we demonstrate that intergenic regulatory regions are an important source of nearly neutral mutations. These results improve our understanding of selection on noncoding regions and the role of nearly neutral mutations for evolutionary processes in outcrossing Brassicaceae species.

Genomic basis of parallel adaptation varies with divergence in Arabidopsis and its relatives.

Parallel adaptation provides valuable insight into the predictability of evolutionary change through replicated natural experiments. A steadily increasing number of studies have demonstrated genomic parallelism, yet the magnitude of this parallelism varies depending on whether populations, species, or genera are compared. This led us to hypothesize that the magnitude of genomic parallelism scales with genetic divergence between lineages, but whether this is the case and the underlying evolutionary processes remain unknown. Here, we resequenced seven parallel lineages of two Arabidopsis species, which repeatedly adapted to challenging alpine environments. By combining genome-wide divergence scans with model-based approaches, we detected a suite of 151 genes that show parallel signatures of positive selection associated with alpine colonization, involved in response to cold, high radiation, short season, herbivores, and pathogens. We complemented these parallel candidates with published gene lists from five additional alpine Brassicaceae and tested our hypothesis on a broad scale spanning ∼0.02 to 18 My of divergence. Indeed, we found quantitatively variable genomic parallelism whose extent significantly decreased with increasing divergence between the compared lineages. We further modeled parallel evolution over the Arabidopsis candidate genes and showed that a decreasing probability of repeated selection on the same standing or introgressed alleles drives the observed pattern of divergence-dependent parallelism. We therefore conclude that genetic divergence between populations, species, and genera, affecting the pool of shared variants, is an important factor in the predictability of genome evolution.

Unboxing mutations: Connecting mutation types with evolutionary consequences.

A key step in understanding the genetic basis of different evolutionary outcomes (e.g., adaptation) is to determine the roles played by different mutation types (e.g., SNPs, translocations and inversions). To do this we must simultaneously consider different mutation types in an evolutionary framework. Here, we propose a research framework that directly utilizes the most important characteristics of mutations, their population genetic effects, to determine their relative evolutionary significance in a given scenario. We review known population genetic effects of different mutation types and show how these may be connected to different evolutionary outcomes. We provide examples of how to implement this framework and pinpoint areas where more data, theory and synthesis are needed. Linking experimental and theoretical approaches to examine different mutation types simultaneously is a critical step towards understanding their evolutionary significance.

On the origin of the widespread self-compatible allotetraploid Capsella bursa-pastoris (Brassicaceae).

Polyploidy, or whole-genome duplication, is a common speciation mechanism in plants. An important barrier to polyploid establishment is a lack of compatible mates. Because self-compatibility alleviates this problem, it has long been hypothesized that there should be an association between polyploidy and self-compatibility (SC), but empirical support for this prediction is mixed. Here, we investigate whether the molecular makeup of the Brassicaceae self-incompatibility (SI) system, and specifically dominance relationships among S-haplotypes mediated by small RNAs, could facilitate loss of SI in allopolyploid crucifers. We focus on the allotetraploid species Capsella bursa-pastoris, which formed ~300 kya by hybridization and whole-genome duplication involving progenitors from the lineages of Capsella orientalis and Capsella grandiflora. We conduct targeted long-read sequencing to assemble and analyze eight full-length S-locus haplotypes, representing both homeologous subgenomes of C. bursa-pastoris. We further analyze small RNA (sRNA) sequencing data from flower buds to identify candidate dominance modifiers. We find that C. orientalis-derived S-haplotypes of C. bursa-pastoris harbor truncated versions of the male SI specificity gene SCR and express a conserved sRNA-based candidate dominance modifier with a target in the C. grandiflora-derived S-haplotype. These results suggest that pollen-level dominance may have facilitated loss of SI in C. bursa-pastoris. Finally, we demonstrate that spontaneous somatic tetraploidization after a wide cross between C. orientalis and C. grandiflora can result in production of self-compatible tetraploid offspring. We discuss the implications of this finding on the mode of formation of this widespread weed.

The Genomic Architecture and Evolutionary Fates of Supergenes.

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

Hybrid seed incompatibility in Capsella is connected to chromatin condensation defects in the endosperm.

Hybridization of closely related plant species is frequently connected to endosperm arrest and seed failure, for reasons that remain to be identified. In this study, we investigated the molecular events accompanying seed failure in hybrids of the closely related species pair Capsella rubella and C. grandiflora. Mapping of QTL for the underlying cause of hybrid incompatibility in Capsella identified three QTL that were close to pericentromeric regions. We investigated whether there are specific changes in heterochromatin associated with interspecific hybridizations and found a strong reduction of chromatin condensation in the endosperm, connected with a strong loss of CHG and CHH methylation and random loss of a single chromosome. Consistent with reduced DNA methylation in the hybrid endosperm, we found a disproportionate deregulation of genes located close to pericentromeric regions, suggesting that reduced DNA methylation allows access of transcription factors to targets located in heterochromatic regions. Since the identified QTL were also associated with pericentromeric regions, we propose that relaxation of heterochromatin in response to interspecies hybridization exposes and activates loci leading to hybrid seed failure.

The genome of Draba nivalis shows signatures of adaptation to the extreme environmental stresses of the Arctic.

The Arctic is one of the most extreme terrestrial environments on the planet. Here, we present the first chromosome-scale genome assembly of a plant adapted to the high Arctic, Draba nivalis (Brassicaceae), an attractive model species for studying plant adaptation to the stresses imposed by this harsh environment. We used an iterative scaffolding strategy with data from short-reads, single-molecule long reads, proximity ligation data, and a genetic map to produce a 302 Mb assembly that is highly contiguous with 91.6% assembled into eight chromosomes (the base chromosome number). To identify candidate genes and gene families that may have facilitated adaptation to Arctic environmental stresses, we performed comparative genomic analyses with nine non-Arctic Brassicaceae species. We show that the D. nivalis genome contains expanded suites of genes associated with drought and cold stress (e.g., related to the maintenance of oxidation-reduction homeostasis, meiosis, and signaling pathways). The expansions of gene families associated with these functions appear to be driven in part by the activity of transposable elements. Tests of positive selection identify suites of candidate genes associated with meiosis and photoperiodism, as well as cold, drought, and oxidative stress responses. Our results reveal a multifaceted landscape of stress adaptation in the D. nivalis genome, offering avenues for the continued development of this species as an Arctic model plant.

HEARTBREAK Controls Post-translational Modification of INDEHISCENT to Regulate Fruit Morphology in Capsella.

Morphological variation is the basis of natural diversity and adaptation. For example, angiosperms (flowering plants) evolved during the Cretaceous period more than 100 mya and quickly colonized terrestrial habitats [1]. A major reason for their astonishing success was the formation of fruits, which exist in a myriad of different shapes and sizes [2]. Evolution of organ shape is fueled by variation in expression patterns of regulatory genes causing changes in anisotropic cell expansion and division patterns [3-5]. However, the molecular mechanisms that alter the polarity of growth to generate novel shapes are largely unknown. The heart-shaped fruits produced by members of the Capsella genus comprise an anatomical novelty, making it particularly well suited for studies on morphological diversification [6-8]. Here, we show that post-translational modification of regulatory proteins provides a critical step in organ-shape formation. Our data reveal that the SUMO protease, HEARTBREAK (HTB), from Capsella rubella controls the activity of the key regulator of fruit development, INDEHISCENT (CrIND in C. rubella), via de-SUMOylation. This post-translational modification initiates a transduction pathway required to ensure precisely localized auxin biosynthesis, thereby facilitating anisotropic cell expansion to ultimately form the heart-shaped Capsella fruit. Therefore, although variation in the expression of key regulatory genes is known to be a primary driver in morphological evolution, our work demonstrates how other processes-such as post-translational modification of one such regulator-affects organ morphology.

Differential Expression of Immune Genes between Two Closely Related Beetle Species with Different Immunocompetence following Attack by Asecodes parviclava.

Endoparasitoid wasps are important natural enemies of many insect species and are major selective forces on the host immune system. Despite increased interest in insect antiparasitoid immunity, there is sparse information on the evolutionary dynamics of biological pathways and gene regulation involved in host immune defense outside Drosophila species. We de novo assembled transcriptomes from two beetle species and used time-course differential expression analysis to investigate gene expression differences in closely related species Galerucella pusilla and G. calmariensis that are, respectively, resistant and susceptible against parasitoid infection by Asecodes parviclava parasitoids. Approximately 271 million and 224 million paired-ended reads were assembled and filtered to form 52,563 and 59,781 transcripts for G. pusilla and G. calmariensis, respectively. In the whole-transcriptome level, an enrichment of functional categories related to energy production, biosynthetic process, and metabolic process was exhibited in both species. The main difference between species appears to be immune response and wound healing process mounted by G. pusilla larvae. Using reciprocal BLAST against the Drosophila melanogaster proteome, 120 and 121 immune-related genes were identified in G. pusilla and G. calmariensis, respectively. More immune genes were differentially expressed in G. pusilla than in G. calmariensis, in particular genes involved in signaling, hematopoiesis, and melanization. In contrast, only one gene was differentially expressed in G. calmariensis. Our study characterizes important genes and pathways involved in different immune functions after parasitoid infection and supports the role of signaling and hematopoiesis genes as key players in host immunity in Galerucella against parasitoid wasps.

Population Genomics of Transitions to Selfing in Brassicaceae Model Systems.

Many plants harbor complex mechanisms that promote outcrossing and efficient pollen transfer. These include floral adaptations as well as genetic mechanisms, such as molecular self-incompatibility (SI) systems. The maintenance of such systems over long evolutionary timescales suggests that outcrossing is favorable over a broad range of conditions. Conversely, SI has repeatedly been lost, often in association with transitions to self-fertilization (selfing). This transition is favored when the short-term advantages of selfing outweigh the costs, primarily inbreeding depression. The transition to selfing is expected to have major effects on population genetic variation and adaptive potential, as well as on genome evolution. In the Brassicaceae, many studies on the population genetic, gene regulatory, and genomic effects of selfing have centered on the model plant Arabidopsis thaliana and the crucifer genus Capsella. The accumulation of population genomics datasets have allowed detailed investigation of where, when and how the transition to selfing occurred. Future studies will take advantage of the development of population genetics theory on the impact of selfing, especially regarding positive selection. Furthermore, investigation of systems including recent transitions to selfing, mixed mating populations and/or multiple independent replicates of the same transition will facilitate dissecting the effects of mating system variation from processes driven by demography.

Impact of demography on linked selection in two outcrossing Brassicaceae species.

Genetic diversity is shaped by mutation, genetic drift, gene flow, recombination, and selection. The dynamics and interactions of these forces shape genetic diversity across different parts of the genome, between populations and species. Here, we have studied the effects of linked selection on nucleotide diversity in outcrossing populations of two Brassicaceae species, Arabidopsis lyrata and Capsella grandiflora, with contrasting demographic history. In agreement with previous estimates, we found evidence for a modest population size expansion thousands of generations ago, as well as efficient purifying selection in C. grandiflora. In contrast, the A. lyrata population exhibited evidence for very recent strong population size decline and weaker efficacy of purifying selection. Using multiple regression analyses with recombination rate and other genomic covariates as explanatory variables, we can explain 47% of the variance in neutral diversity in the C. grandiflora population, while in the A. lyrata population, only 11% of the variance was explained by the model. Recombination rate had a significant positive effect on neutral diversity in both species, suggesting that selection at linked sites has an effect on patterns of neutral variation. In line with this finding, we also found reduced neutral diversity in the vicinity of genes in the C. grandiflora population. However, in A. lyrata no such reduction in diversity was evident, a finding that is consistent with expectations of the impact of a recent bottleneck on patterns of neutral diversity near genes. This study thus empirically demonstrates how differences in demographic history modulate the impact of selection at linked sites in natural populations.

Genetic basis and timing of a major mating system shift in Capsella.

A crucial step in the transition from outcrossing to self-fertilization is the loss of genetic self-incompatibility (SI). In the Brassicaceae, SI involves the interaction of female and male specificity components, encoded by the genes SRK and SCR at the self-incompatibility locus (S-locus). Theory predicts that S-linked mutations, and especially dominant mutations in SCR, are likely to contribute to loss of SI. However, few studies have investigated the contribution of dominant mutations to loss of SI in wild plant species. Here, we investigate the genetic basis of loss of SI in the self-fertilizing crucifer species Capsella orientalis, by combining genetic mapping, long-read sequencing of complete S-haplotypes, gene expression analyses and controlled crosses. We show that loss of SI in C. orientalis occurred < 2.6 Mya and maps as a dominant trait to the S-locus. We identify a fixed frameshift deletion in the male specificity gene SCR and confirm loss of male SI specificity. We further identify an S-linked small RNA that is predicted to cause dominance of self-compatibility. Our results agree with predictions on the contribution of dominant S-linked mutations to loss of SI, and thus provide new insights into the molecular basis of mating system transitions.

Long-term balancing selection drives evolution of immunity genes in Capsella.

Genetic drift is expected to remove polymorphism from populations over long periods of time, with the rate of polymorphism loss being accelerated when species experience strong reductions in population size. Adaptive forces that maintain genetic variation in populations, or balancing selection, might counteract this process. To understand the extent to which natural selection can drive the retention of genetic diversity, we document genomic variability after two parallel species-wide bottlenecks in the genus Capsella. We find that ancestral variation preferentially persists at immunity related loci, and that the same collection of alleles has been maintained in different lineages that have been separated for several million years. By reconstructing the evolution of the disease-related locus MLO2b, we find that divergence between ancient haplotypes can be obscured by referenced based re-sequencing methods, and that trans-specific alleles can encode substantially diverged protein sequences. Our data point to long-term balancing selection as an important factor shaping the genetics of immune systems in plants and as the predominant driver of genomic variability after a population bottleneck.